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Уточнена математическая модель, описывающая процессы переноса и ассимиляции вредных веществ в по- 
мещениях участков обкатки в отношении выбросов твердых частиц: сажи и бенз(а)пирена. Модель позволя- 
ет прогнозировать качество воздуха в рабочей зоне таких помещений, определять опасные канцерогенные 
зоны в поле концентраций токсичных газов. 
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Введение. Авторами была разработана нестационарная трехмерная математическая модель кон- 
вективно-диффузионного распространения вредных веществ в производственном помещении, 
которая позволила достаточно точно определить поля концентраций в отношении теплых газооб- 
разных токсичных веществ, являющихся легкими примесями (аэрозолями) [1]. Однако вопрос 
распространения и ассимиляции твердых частиц, а именно частиц сажи и бенз(а)пирена, присут- 
ствующих в отработанных газах (ОГ), отходящих от стендов обкатки комбайнов, не рассматривался. 
Постановка задачи. Известно, что канцерогенные вещества являются наиболее токсичными 
компонентами ОГ. Ряд коэффициентов относительной агрессивности токсичных компонентов ОГ 
приведен в работе [2]: 

Ксо : Ксн : Кмох : Ктч : Кзох : КБал = 1: 3,16 : 41,1:41,1 : 16,5 : 1260000. (1) 

По сравнению с угарным газом бенз(а)пирен опаснее для человека более чем в миллион 
раз. Из дисперсионного анализа известно, что дизельная сажа по массе представлена твердыми 
частицами размера от 0,7 до 15 мкм, которые наиболее опасны для здоровья человека, так как 
попадают через органы дыхания и переносят частички бенз(а)пирена. 

Поэтому расчетные опасные зоны по канцерогенам могут существенно отличаться от зон 
загазованности, соответственно снижается точность прогнозирования качества воздуха в иссле- 
дуемых производственных помещениях [1]. Поэтому предпринята попытка уточнить опасные зоны 
загазованности, рассмотрев распространение тяжелых примесей. 

Для решения поставленной задачи необходимо уточнение ранее разработанной матема- 
тической модели [1] как в отношении интенсивности выделения бенз(а)пирена источником за- 
грязнения (двигатель машины), так и в отношении особенностей распространения твердых час- 
тиц, которые определяются скоростью их осаждения в воздушной среде помещения обкатки. 
Теоретические положения. Математическая модель переноса. Предложенный метод рас- 
чета базируется на созданной численной модели [1] процесса миграции примеси в помещении на 
основе трехмерного уравнения конвективно-диффузионного переноса вредностей в воздушной 
среде 

дС/0! + 9(иС) [ах + 9(уС)/ду + 9ВС(и- м.) /62 = 


= 9/0х(А, ЭС/ 9х) + 9/ду(А, ОС/бу) + 0/02 (А, ОС] 2) + У` 0, (06 — ›), 


где С - концентрация загрязняющего вещества в помещении; и‚у,и»— компоненты вектора 
скорости воздушной среды; и. — скорость подъема/оседания легких/тяжелых примесей; 


5 
А - коэффициент турбулентного обмена;О;(1) — интенсивность выброса вредного вещества 
от 7-го источника в помещении; д0(/—/;) - дельта-функция Дирака; /; = (х, у, 2) - координаты 
источника выброса. 
Это уравнение неразрывности потока примеси. Члены, содержащие компоненты скорости 
и,у,и, в направлении координатных осей х, у, 7, описывают конвективный перенос примеси. 
В правой части уравнения сгруппированы члены, описывающие диффузионную часть. 


(2) 


534 


Вестник ДГТУ, 2010. Т. 10. №4(47) 


Для расчета поля скорости воздушного потока в помещении, индуцированного работой 
вентиляции, делается допущение, что движение воздушной среды в помещении потенциальное, 
тогда компоненты скорости воздушной среды определяются соотношениями: 

и=дР/ох, у=дР/ду, и=дР/02, 
где Р- потенциал. 
Уравнение для определения потенциала имеет вид 
0?Р/дх? +0?Р/ду? +6?Р/62? =0. (3) 

Для решения дифференциальных уравнений (2) и (3) задаются начальные и граничные 
условия [1]. 

Распространение тяжелой примеси. Особенности распространения тяжелых примесей опре- 
деляются в значительной мере собственной скоростью их осаждения, которая зависит от плотно- 
сти, размеров и формы частиц аэрозоля. 

Скорость падения/оседания сферических твердых частиц может быть определена по фор- 
муле Стокса, полученной из условия равновесия, действующего на шаровую частицу трех сил 
(тяжести, сопротивления и поддерживающей): 

и, =1,3.10-2р‚2, (4) 
где р; — плотность частиц тяжелой примеси, //— радиус частицы. 

Таким образом, для определения скорости оседания твердой частицы примеси необходи- 

мо знание физических свойств вещества и данных дисперсионного анализа. 
Свойства и дисперсионный анализ дизельной сажи. Изучением дисперсионного состава 
дизельной сажи занимались многие исследователи [3, 4]. В их работах можно найти подробные 
ссылки на современные исследования в данной области. Утверждения авторов этих исследований 
зачастую полярно противоположные. В некоторых указывается, что средний размер частиц в вы- 
хлопных газах дизельного двигателя колеблется в пределах 0,2-0,3 мкм и слабо зависит от режи- 
ма работы и типа исследуемого двигателя [3]. В других приведены распределения твердых частиц 
в выхлопных газах в зависимости от режима работы двигателя и отмечено, что размеры частиц 
меняются в пределах 1-14 мкм [4]. 

Столь разные результаты не могут быть отнесены к погрешностям измерений и нет при- 
чины сомневаться в корректности проведения экспериментов различными исследователями. Дело 
в том, что классический процесс конденсации начинается с флуктуаций плотности газового объе- 
ма с последующим ростом образовавшихся микрочастиц [5]. Поэтому, строго говоря, в распреде- 
лении частиц по размерам должны присутствовать частицы от молекулярных до макроскопиче- 
ских размеров. Размер микрочастиц определяется конкретными экспериментальными условиями, 
а универсальной причины ограничения роста частиц просто не существует. На основе разрабо- 
танного метода голографической регистрации конденсированной фазы в ОГ дизельного 
двигателя проведены измерения характеристик дисперсного состава ОГ дизельного двигателя 
КамАЗ-740 при различных режимах работы и изменения дисперсного состава (рис.1). 
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Рис.1. Распределение массы (объема) частиц сажи 
по диаметрам до каталитического нейтрализатора 
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Особенностью метода является использование для регистрации голограммы скоростной 
цифровой камеры с прямым вводом в ПЭВМ и непосредственной компьютерной обработки полу- 
чаемых результатов. Как видно из гистограммы, наибольшую массу имеют частицы двух харак- 
терных диаметров: 5 мкм и 15 мкм. Процесс их распространения можно описывать на основе 
уравнения (2) модели массопереноса, каждый раз уточняя величину и’ (4). При этом воспользу- 
емся упрощением, что вся масса выбрасываемых частиц имеет диаметр 10 мкм. 

Характеристика интенсивности движения газовоздушных потоков и примесей в про- 
изводственном помещении. Большое влияние на поля концентраций вредные вещества (ВВ) 
оказывают коэффициенты турбулентной диффузии, входящие в уравнение массопереноса (2). 

При турбулентном течении потока и турбулентном переносе коэффициенты Ал пропорцио- 
нальны диссипируемой энергии и определяющему размеру рассматриваемого источника. Для вен- 
тилируемых помещений, зная величину энергии, поступающей в помещение, и определяющий 
размер источника, можно определить коэффициент турбулентного обмена в плоскости (х, у), м*/с: 


и 4 
ЮО и, (5) 
где & —- кинетическая энергия воздушных потоков в помещении, 1 - определяющий размер поме- 


щения для прямоугольных отверстий и сечений неправильной формы, /[„ = (ЁЕ- пло- 


щадь сечения помещения, перпендикулярная направлению движения воздушных потоков). 
Коэффициент турбулентного обмена по высоте помещения находится из зависимости, 
предложенной М.Е. Берляндом для наибольших классов устойчивости воздушной среды 


Апр = А; _ ан что допустимо для исследуемого помещения. Здесь А?-, —- коэффициент турбу- 


лентного обмена на высоте 1 м; 2-— текущая высота помещения, м [6]. 
Кинетическая энергия воздушных потоков в помещении определяется из формулы: 


где пс - энергия приточных струй; = -— тепловых струй; =п - движущихся предметов; 

&вс — вытяжных струй. 

Величины &цс И &тс учитываются в модели рассеивания ВВ: 
= И 

енс = СъёУ? /2С , (6) 
где Св- масса воздуха, подаваемого в помещение, кг/с, Сь =[р, где [ — объемный расход воздуха, 
м^/с; р — плотность воздуха, кг/м?; & —- коэффициент местного сопротивления на выходе из 
воздухораспределителей (табличные справочные данные); и -— средняя скорость выхода 


воздуха из приточных отверстий, м/с; Сл -— масса воздуха в объеме помещения, кг, (б- =\ф). 
Учитывая, что кратность воздухообмена к, = Г/И , где Г’ - свободный объем помеще- 


ния, м?, после преобразований получаем: 





а а №. (7) 
р —- ее 

Стс = тс _ 802 1+п ИЛИ Стс — К 2, (8) 
Сп С „ор ь Рв 


где О - теплонапряженность объема, Вт/м?, О=Оз/Ит; Оз - теплопоступление явного тепла в 


помещение, Вт; © — ускорение свободного падения, м“/с; 2 — расстояние от полюса до рас- 
сматриваемой точки, м: 2 =2, + 2и, ГДе 2, - расстояние от источника до полюса струи, м; 
2„ — расстояние от источника до рассматриваемой точки, м; п — экспериментальная констан- 
та, п=0,5; С› - теплоемкость воздуха в удалении от источника; То — температура воздуха в 
удалении от источника, К; р, — плотность воздуха в удалении от источника, кг/м; А — коэф- 
фициент пропорциональности для нагретой пластины, расположенной на основании, 
к= 3,43. 10°, расположенной заподлицо с плоскостью, А=18,3. 10°. 

Явное теплопоступление в помещение учитывается для источников, оборудованных мест- 

ными отсосами: 


Оя = СогОог (Тог -Т, И -Фог) Ад, 
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а для источников, ими не оборудованных: 
хх эхх ип-ХХ 
О; = СрогОог (Гог -Т5)-Фог) Ад, 
где Сьог, о. — изобарные теплоемкости ОГ на режимах максимальной нагрузки и холостого хо- 


тах 


да, Дж/кг-К; Того, р. — температура ОГ на различных режимах нагрузки, К; Ад — количе- 


ство стендов; ф — степень эффективности работы местного отсоса. 
Характеристика источника выделения вредностей. Когда двигатель комбайна работает под 
нагрузкой, с ОГ через неплотность местных отсосов в штатных условиях или при разрыве отсоса в 
аварийной ситуации в воздух РЗ поступают токсичные компоненты. Для определения их массово- 
го выброса была разработана эмпирическая модель на основе регрессионных полиномов [7]. 
Функция ОЕ) может быть определена по формуле: 


(0 =Сл0 _-Фог) /3600, (9) 


где с. - концентрация ВВ в ОГ от #го источника загрязнения, г/м; Оотерд — расход отрабо- 


Д 
танных газов двигателя комбайна при средней мощности обкатки, м?/с; О эффектив- 


ность местного отсоса у обкаточного стенда, которая либо рассчитывается, либо принима- 
ется равной от 70 до 90 % в зависимости от его конструктивного устройства. 

Концентрации ВВ в ОГ для сажи и бенз(а)пирена могут быть вычислены по полиномам, 
полученным в результате аппроксимации большого числа достоверных экспериментальных дан- 
ных по выбросам этих веществ (табл.1). 

Таблица 1 
Концентрации вредных веществ в ОГ комбайновых ДВС в зависимости 
от режима нагрузки (установившиеся режимы работы двигателя при обкатке) 


Концентрации ВВ в ОГ источника загрязнения в зависимости от его относительной мощности Точность расчетов, А 
Сс=0,2114№2 —0,0214№+0,0324 , г/м? 0,9838 


Сон, =-2,2933№3 +5,3143 №? +-0,2290№+0,4313 МКт/м: 





Для приближенных оценочных расчетов можно воспользоваться данными удельных пока- 
зателей выбросов сажи и бенз(а)пирена дизельными двигателями при обкатке из нормативной 


литературы. 
Расход ОГ на любом режиме обкатки комбайна (м?/с), определяется по формуле: 
Оо. = 24,353ф „В Ца М, [НПт„ р. , (10) 
где ра коэффициент продувки, для 2-тактных двигателей Ф:0 =1,4-1,6; для 4-тактных 


Ф.0 =1,1-1,25; В- коэффициент молекулярного изменения рабочей смеси, определяемый из 


теплового расчета; а — коэффициент избытка воздуха; Н,, — низшая теплота сгорания топ- 
лива, Дж/кг; №, — эфективная мощность двигателя, кВт; Г, — теоретический расход воздуха 
на сгорание, кмоль; и; 7; — индикаторный и механический КПД двигателя. 

Плотность смеси ОГ вычисляли для смеси осредненного состава на режимах прогрева, хо- 
лостом ходу и при пробеге. 

Численное решение уравнений (2) и (3) осуществляется с помощью метода конечных эле- 
ментов. При решении уравнений предложенной математической модели использовались процеду- 
ры из пакета НехРОЕ — программы, поддерживающей метод конечных элементов при моделиро- 
вании объектов с распределенными переменными, описываемыми нелинейными дифференциаль- 
ными уравнениями с частными производными. 

Результаты модельных расчетов. Объект исследования - участок обкатки комбайнов сбороч- 
ного цеха № 20 ОАО «КЗ “Ростсельмаш”». Исходные данные для модельных расчетов: геометри- 
ческие параметры помещения обкатки комбайнов (длина - 132 м, ширина - 60 м, высота - 14 м); 
количество и геометрические размеры стендового оборудования и обкатываемых сборочных еди- 
ниц; время работы в соответствии с технологическим процессом обкатки; виды и расположение 
источников загрязнения (обкатываемых машин) и рабочих мест; степень эффективности местных 
отсосов (70 %); интенсивность движения газовоздушных потоков (входные и выходные скорости); 
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интенсивность выделения ВВ от источников, полученных по регрессионным полиномам в зависи- 
мости от режима обкатки (табл.2); схема организации движения газовоздушных потоков в венти- 
лируемом помещении с различным расположением воздухораспределителей системы механиче- 
ской вентиляции, а также приточных и вытяжных отверстий (рис.2); расчетный период года - те- 
ПЛЫЙ И ХОЛОДНЫЙ. 
Таблица 2 
Интенсивность выделения сажи и бенз(а)пирена источником загрязнения 


Максимально разовые выбросы вредных веществ при режимах обкатки 
комбайна АСКО5-530 
о Бенз(а)пирен, г/с 
обкатки , , 


Штатный Авария Неисправность Штатный Авария Неисправность 
режим, (отрыв обкатываемой режим, (отрыв обкатываемой 


машины 


ф=0,7 отсоса) машины ф=0,7 отсоса) машины 


0,08706 0,2902 8,01 : 107 2,67: 10° 8,01 : 10° 


т 0,0717 0,2151 0,15: 10° 0,5- 10° р 


нагрузки 


[Холостой ход 1 00610" 0,18 10° 


выход 











Вариант 1 


выходные 
отверстия 


(вентиляторы) 


ты комбайны 
© 
Рис.2. Расчетная схема организации движения газовоздушных потоков 
в помещении обкатки комбайнов 


При проведении модельных расчетов были рассмотрены как штатные режимы работы 
оборудования на участке обкатки комбайнов, так и различные аварийные ситуации, которые ока- 
зывают существенное влияние на формирование полей концентрации ВВ и опасных зон загазо- 
ванности в помещении. 

Под штатным режимом понимается исправная работа всех видов оборудования и систем 
вентиляции производственного помещения. Загазованность воздуха рабочей зоны создается лишь 
за счет выделения ВВ от источников, оборудованных негерметичными местными отсосами, а так- 
же мест проверки работы кондиционеров комбайнов, которые не оснащены системой локальной 
вытяжной вентиляции. Таким образом, загазованность определяется фоном, полученным от всех 
источников загрязнения. 

Под аварийными ситуациями на участке обкатки машин понимаются вероятные случаи 
разрыва местных отсосов, которыми оборудованы все стенды обкатки комбайнов, а также ситуа- 
ция, когда внезапно выключается общеобменная и аварийная механическая вентиляция. Это наи- 
более типичные аварийные ситуации для данного производства. Рассмотрены наиболее вероят- 
ные и неблагоприятные случаи, причем моделировалось и срабатывание систем защиты (ручное и 
автоматическое отключение оборудования, включение аварийной вентиляции и т.д.). 

Штатный режим работы участка обкатки. Численные эксперименты на модели. Расчет- 
ная схема организации движения газовоздушных потоков в помещении обкатки комбайнов соот- 
ветствует «Варианту 1» (реализован в ОАО «КЗ “Ростсельмаш”») (рис.2). Исправно функциониру- 
ют системы принудительной общеобменной и местной вентиляций помещения. Поля концентра- 
ций бенз(а)пирена в воздухе помещения участка обкатки комбайнов получены при нормальной 
работе местных отсосов (фот = 0,7) для двух расчетных периодов года (теплого и холодного). 
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Значения входной и выходной скорости движения воздуха в вентилируемом помещении принима- 
лись в диапазоне, установленном нормативными значениями микроклиматических параметров — 
0,2 м/с. Поля концентраций бенз(а)пирена рассчитаны по модели на трех характерных срезах по 
высоте помещения 2, 4, 12 м, что обусловлено реальным расположением разных рабочих мест на 
участке. В рассматриваемом случае, загазованность производственного помещения создается за 
счет выделения газов и бенз(а)пирена от семи источников загрязнения, которые определяют фо- 
новые концентрации в воздухе, увеличиваясь с течением времени. 





Е 

90 
1.15 
1.10 
1.05 
1.00 
0.95 
60 0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
30 0.60 
0.55 
0.50 
0.45 
0.40 
0.35 
0.30 
о 0.25 
0.20 
0.15 
0.10 
0.05 
— 0.00 

—30` 





[_] - СОНИТОВНО-2И2ИЕНИЧеСКаЯ 


0 30 60 90 120 Х 
зона (по бенз@/пирениу/ 


а) 





Сопсепегайоп 











® 


6) 


Рис.3. Расчет значений концентраций бенз(а)пирена в воздухе участка обкатки в холодный период года: 
а — поля относительных значений концентраций бенз(а)пирена; 
6 — динамика изменения относительных концентраций бенз(а)пирена 


На рис.3,а представлены рассчитанные по модели поля относительных концентраций 
бенз(а)пирена на высоте 4м на один час работы участка в холодный расчетный период года. Мак- 
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симальная концентрация бенз(а)пирена - х. Как следует из картин распределения полей концен- 
трации, существенное влияние на их формирование оказывает расчетный период года. Объяс- 
няется это влияние тем, что в холодный период года в соответствии с технологией обкатки в те- 
чение рабочего часа выездные ворота остаются закрытыми. Результаты модельных расчетов по- 
лей подвижности воздуха в помещении участка обкатки подтверждают, что в холодный период 
года скорости движения газовоздушных потоков в правой части помещения (у выездных ворот) 
невелики, что создает условия для образования застойной зоны и, как следствие, повышения 
концентрации в этой части помещения (рис.3,а). Как видно из результатов модельных расчетов, 
концентрация по бенз(а)пирену превышает установленные предельно допустимые концентрации 
(ПДК) на 15 % даже в штатном безаварийном режиме работы обкаточного цеха. При открытых 
въездных и выездных воротах в теплый период года создают хорошие условия для дополнитель- 
ного проветривания, и превышения нет. 

Авария, связанная с вероятным выходом из строя общеобменной вентиляции. Рассмот- 
рим наиболее нежелательную аварийную ситуацию, когда внезапно отключается (не работает) 
принудительная общеобменная и аварийная вентиляция. Модельные расчеты проведены для теп- 
лого и холодного периода года при существующем на предприятии варианте организации возду- 
хообмена помещения. 

На рис.4,а, 6 показаны поля относительной концентрации по бенз(а)пирену. Расчеты про- 
изведены при скорости движения воздуха в помещении 0,2 м/с и эффективности местных отсосов 
70 %. Следует также отметить, что поля сформированы на час работы, т.е. на максимальное вре- 
мя продолжительности аварийной ситуации. 
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Рис.4. Поля относительных значений концентраций бенз(а)пирена в воздухе участка обкатки 
в аварийном режиме: а - теплый период года; 6 — холодный период года 
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Анализ расчетных полей концентрации (рис.4,а, 6) показывает, что при реализации схемы 
движения газовоздушных потоков по «Варианту 1» опасная канцерогенная зона по бенз(а)пирену 
в холодный период года достигает почти четверти от площади всего помещения и представляет 
собой значительную опасность для работающих. В этой канцерогенной зоне три рабочих места, 
а ПДК превышен более чем в два раза. В теплый период года канцерогенные зоны имеют локаль- 
ный характер и находятся около двух рабочих мест, но санитарно-гигиенический норматив также 
превышается в 1,5 раза. 

Результаты расчетов доказывают, что в случае неработающей обменной и аварийной вен- 

тиляции (а для ряда предприятий России это, к сожалению, штатная ситуация) операторы обка- 
точных участков находятся в экстремально-негативных условиях. 
Авария, связанная с вероятной разгерметизацией местного отсоса. Рассмотрим внезап- 
ный «разрыв» крайнего местного отсоса на стенде. На рис.5 показаны картины полей концентра- 
ций бенз(а)пирена в воздухе участка обкатки, которые сформировались за разное время. Началь- 
ный период работы был штатным (около 40 мин), затем — фаза периода аварии та. 
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Рис.5. Расчет значений концентраций бенз(а)пирена в воздухе участка обкатки в холодный период года: 
а — динамика изменения в течение всего аварийного периода, разрыв крайнего источника, ручное отключение; 
6 — поля относительных значений концентрации бенз(а)пирена в воздухе участка обкатки во внештатной ситуации 
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Как видно, при ручном отключении поврежденного стенда и оборудования т. = 300 с 
кратковременно формируется локальная канцерогенная зона бенз(а)пирена (рис.5). После отклю- 
чения оборудования превышение ПДК концентрации не наблюдается. 

Нештатная ситуация, вызванная неисправностью обкатываемых комбайнов. На участке 
обкатки имеют место случаи, когда один или несколько комбайнов неисправны или не отрегули- 
рованы. Выбросы вредных веществ таких комбайнов превышают «паспортные» значения в десять 
и более раз. Например, комбайны 2 и 5 в холодный период года неисправны (превышение значе- 
ний по бенз(а)пирену в 10 раз). За один час работы комбайнов площадь, где имеет место пре- 
вышение ПДК, составила 88,63 % от всей площади обкаточного участка (рис.5,6). 

Результаты расчетов и их обсуждение. Полученные результаты при дальнейших исследова- 
ниях могут быть уточнены в отношении: 

— учета вихревых движений воздушных токов в помещении при расчетах активной венти- 
ляции в загазованных помещениях; 

— применения дискретного подхода к расчету распространения твердых частиц с учетом 
их формы, размера и траектории движения. 

Также можно получить поля распределения концентраций сажи в воздухе рабочей зоны 
в объеме всего помещения. 

Выводы. Уточнена математическая модель процессов переноса и ассимиляции вредных веществ 
в воздушной среде участка обкатки комбайнов в отношении тяжелых примесей. 

Получены поля концентрации бенз(а)пирена в воздухе помещения участка обкатки ком- 
байнов, которые позволили определить канцерогенные зоны по бенз(а)пирену, негативно воздей- 
ствующие на работников в условиях активной вентиляции. 
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